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METEOROLOGICAL ASPECTS OF STORM SURGE GENERATION?@ 


D. Lee Harris! 
(Proc. Paper 1859) 


INTRODUCTION 


The coastal floods, produced by Hurricane AUDREY in June 1957 claimed 
more than 300(1) lives and destroyed more than 100 homes. Even in this en- 
lightened year of 1957, when hurricane warnings are broadcast by radio and 
television long before the storm arrives, the high tides produced by the storm 
continue to claim more lives than all other aspects of hurricane behavior. Al- 
though many of the newcomers in the Cameron area evacuated the threatened 
region soon after the hurricane warnings were announced, many of the old- 
timers remained because they had weathered previous hurricanes with little 
difficulty and they did not expect this one to be more severe than the others. 

A few facts concerning hurricanes and their effects on sea level, which if fully 
understood, may prevent future disasters of this kind, are presented in this 
paper. 

This paper is concerned primarily with the meteorological aspects of 
storm surge generation by hurricanes. The physical aspects of hurricane 
storm surge behavior along the Atlantic and Gulf coasts of the United States 
have been discussed in greater detail in Harris.(2) Zetler(3) gave an extensive 
treatment of the effects of hurricanes on sea level at Charleston, S. C. 
Redfield and Miller(4) considered the effects of several recent Atlantic Coast 
hurricanes, Hubert and Clark(5) published charts showing the high water 
marks established after 16 Atlantic and Gulf coast hurricanes. Reid(6) de- 


veloped a storm tide index based on the size and intensity of the hurricane and 
the slope of the continental shelf. 


Note: Discussion open until May 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1859 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 7, December, 1958. 


a. Presented at the ASCE Hydraulics Conference in Cambridge, Mass., 
August, 1957. 


Office of Meteorological Research, Storm Surge Research Unit, U. S. 
Weather Bureau, Washington, D. C. 
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The wind and pressure distributions of a hurricane are never known with 
desired accuracy before the storm strikes land. Hence, in forecasting storm 
surges and in the design of engineering works to be constructed as protection 
from some future hurricane, it is necessary to use a storm model. Such 
models can be coustructed in a variety of ways depending on the use for which 
they are intended. For storm surge work it is desired to describe as accu- 
rately as possible the distribution of atmospheric pressure and sustained wind 
speeds near the surface of the sea. 

The pressure in many hurricanes can be represented with a close approxi- 
mation by an equation of the type 


P =P, + (P, - P,) F (r/R), (1) 


where P,, = the pressure just outside the storm, 


Po = the pressure at the center of the storm, 


R =some characteristic radius of the storm, 


r = distance from the center of the storm, 
F(r/R) indicates a function of r/R. 


The radius of the eye, the distance from the center to the zone of maximum 
winds, and the distance from the center to the outermost closed isobar have 
all been suggested for use as the characteristic radius. The proper choice 
probably depends on whether it is desired to obtain the maximum accuracy for 
the high wind speed zone or for the entire storm. (Pp - Po) is intended as a 
measure of the intensity of the storm, but Pp is likewise difficult to define 
explicitly. Some workers prefer to use the pressure of the outermost closed 
isobar, still others prefer some other definition. Virtually all workers in 
this field use the best obtainable estimate of the lowest pressure in the storm 
for Po: Because an objective definition of P, is easily constructed, many 
writers use this value to estimate the intensity of the storm. 

The first approximation of the wind speed within any cyclone is given by 
the gradient wind equation. This can be expressed in terms similar to Eq. (1) 
by: 


Ve /r + tv. 


= P'(P, - (2) 


Vg = gradient wind speed 
f =coriolis parameter 
P = density of air 

F' = derivative of F with respect to (r/R) 


Near the center of a hurricane Vg2/r is considerably larger than fV g? hence 
the second term on the left can be neglected within the accuracy of this” ap- 
proximation. Thus the gradient wind, and presumably, the true surface wind 
also, is approximately proportional to the square root of (P, - Po). This re- 


lationship is often used in estimating the maximum wind speed in a 
hurricane, (11) 
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A hurricane model developed along the lines indicated by Eqs. (1) and (2) is 
described by the following equations: 


Vie + = - (4) 
f 
The symbol R in this model has been identified with the radius of maximum 
wind speed. The agreement between this model and the observations of an ex- 
ceptionally well observed hurricane is shown in Fig. 1. The insert map gives 
the storm track. The dots on the left give the observed pressure at several 
stations in the vicinity of Lake Okeechobee, Florida. The solid line gives the 
theoretical pressure profile fitted to three points within 50 miles of the storm 
center. The corresponding theoretical wind profile is given by the upper curve 
on the right hand side of the figure. The observed winds at one station are 
given by the dots below this curve. A solid line has been drawn through these 
dots by eye in order to obtain a smooth profile. The observed wind speed 
varies in a systematic way from about 65% of the computed wind speed at the 
outer edge of the data shown to almost 90% of the predicted value near the 
zone of maximum wind speed. Equally good agreement between the theoreti- 
cal and observed wind speeds has been obtained in only a few storms. This 
lack of agreement between the theoretical and observed winds is due in part 
to the elementary nature of the model, but perhaps equally to the lack of first 
class wind records near the center of hurricanes. The parameters for this 
model have been computed for a large number of storms and were first 
published in Hydrometeorological Report No. 32.(7) 

Data for new storms are added and data for old storms are corrected as 
new information becomes available. An up-to-date (August 1, 1957) version of 
these data is included in this paper in Table L. 

Eqs. (2) and (4) can be valid only for stationary storms. The wind speeds 
on the right hand side of a moving storm are generally higher than those ob- 
served to the left of the storm track. A simple but incomplete explanation, 
widely accepted for more than half a century, is based on the assumption that 
the forward speed of the storm must be added to the rotary speed on the right 
hand side and subtracted on the left. This would give an asymmetry in the 
wind velocity on the right and left sides of the storm track due to the forward 
speed of the storm. This assumption usually leads to corrections in the right 
direction but not of the proper magnitude. Winds of hurricane intensity are 
significantly affected by friction with the underlying water or land surface. 
The effects of friction and of the inflow necessary to feed the vertical currents, 
prominent throughout the region of hurricane winds, may exceed the effects of 
the forward motion of the storm, at least near the surface. 

Although Eq. (4) is too simple to permit an adequate determination of the 
wind field in a hurricane, it is possible to make use of empirically derived 
relations between the wind speed at nearby locations over water and over land, 
and between the observed wind speed over various types of frictional surfaces 
and the theoretical wind speed to determine a usefully accurate picture of the 
wind field. Figs. 2a and b show the reconstructed pressure and wind fields 
for the New England Hurricane of September 21, 1938.(8) In the derivation of 
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Aug. 28, 1893 
Oct. 1, 1893 

Oct. 13, 1893 
Oct. 2, 1898 

Oct. 31, 1899 
Sept. 8, 1900 
Aug. 14, 1901 
Sept. 11,1903 
June 17, 1906 
Sept.17, 1906 
Sept.27, 1906 
Oct. 18, 1906 
July 21, 1909 
Sept.20, 1909 
Cct. 11, 1909 
Oct. 17, 1910 
Oct. 18, 1910 
Aug. 28, 1911 
Sept. 3, 1913 
Aug. 17, 1915 
Sept.29, 1915 
July 5, 1916 
Aug. 18, 1916 
Oct. 18, 1916 
Sept.28, 1917 
Sept. 9, 1919 


*Sept.14, 1919 
Sept.21, 1920 
June 22, 1921 
Oct. 25, 1921 
Aug. 25, 1924 

®Aug. 26, 1924 
Oct. 19, 1924 

"Oct. 20, 1924 
Dec. 2, 1925 
July 28, 1926 
Aug. 25, 1926 
Sept.18, 1926 
*Sept.20, 1926 
Oct. 20, 1926 
Sept.16, 1928 
June 28, 1929 
Sept.28, 1929 


28.28 
28.22 
28.33 
28.82 
28.70 
27.64 
28.72 
28.84 
28.91 
28.98 
28.50 
28.84 
28.31 
28.94 
28.30 
27.80 
28.33 
28.92 
28.81 
28.01 
27.79 
28.38 
28.00 
28.76 
28.48 
27.44 
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29.61 
29.99 
29.61 
30 .03 
30.49 
29.78 
30.16 
30.12 
29.98 
30.38 
30 .97 
29 .80 
30.27 
30.30 
30.07 
29.19 
29.77 
36.10 
29.98 
29.65 
30.14 
30.03 
30.77 
30.20 
29.88 
29.78 


near 28.0 


28.17 
28.29 
28.70 
28.70 
28.70 
28.83 
28.95 
28.34 
28.31 
27.59 
28.20 
27.52 
27.62 
28.62 
28.15 


29.40 
30.03 
29.59 
30.33 
29.62 
29.82 
29.62 
29.90 
29.91 
30.35 
29.99 
30.13 
29.97 
30.38 
29.97 
30.98 


1.33 
Lave 
1.28 
1.21 
1.79 
2.14 
1.44 
1.28 
1.07 
1.40 
1.57 

1.96 
1.36 
1.77 
1.39 
1 
1.18 
1.17 
1.64 
2.44 
1.65 
2.77 
1.44 
1.40 
2.29 


inches 

-97 
1.86 
1.30 
1.63 

92 

1.57 
2.04 
2.40 
1.93 
2.45 
2.75 
1.35 
1.93 


958 
956 


976 
972 
936 
973 
977 
979 
981 
965 
977 
959 
980 
958 
941 
959 
979 
976 
949 
938 
961 
948 
974 
964 
929 


980 
954 
958 
972 
972 
972 
975 
98 

960 
959 
934 
955 
932 
935 
969 
953 
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1003 45 
1016 60 
1003 44 
1017 41 
1033 61 
1009 73 
1021 48 
1020 43 
1015 36 
1029 42 
1012 53 
1009 32 
1925 66 
1026 45 
1018 60 
989 48 
1008 49 
1019 40 
1015 39 
1004 55 
1021 
1017 56 
1042 94 
1023 
1012 48 
10U? 73 
1013 33 
1017 63 
1002 44 
1027 55 
1003 31 
1010 38 
1003 27 
1013 33 
1013 53 
1028 69 
1016 82 
1020 65 
1015 83 
1029 94 
1015 46 
1019 66 


Table I - Characteristics of United States Hurricanes 
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a 
mb. 
28.28 958 
28.65 970 
28.33 959 
29.12 986 
28.90 979 
28.48 964 
29.42(1)996 
29.47 9238 
29.46 938 
29.50 999 
28.50 965 
29.26 991 
29.00 982 
29.23 990 
23.36 960 
27.80 941 
28.94 980 
29.02 983 
29.36 934 
28.14 953 
28.01 949 
28.38 961 
28.00 948 
28.76 974 
28.51 966 
27.4% 929 
28.65 970 
28.99 982 
29.37 995 
28.29 958 
28.80 975 
28.71 972 
29.10 985 
29.17 988 
28.80 975 
28.31 959 
27.61 935 
28.20 955 
29.16 988 
27.62 935 
29.12 986 
28.18 954 
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Savannah, Ga. 

Moss Point, Miss. 
South Island, S. C. 
Jacksonville, Fla. 
Charleston, S. C. 
Galveston, Tex. 

New Orleans, La. 
Taupa, Fla. 
Jupiter, Fla. 
Columbia, S. C. 
Ship at Scranton, Miss. 
Jupiter, Fla. 

Bey City, Tex. 

New Orleans, La. 
Sand Key, Fla. 


S.S.Jean Nr. Tortugas, Fla. 


Tuanpa, Fla. 
Savannah, Ga. 
Raleigh, N.C. 
Velasco, Tex. 
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New Orleans, La.(Pauline St.Wharr) 12 


Ft. Morgan, Ala. 

Santa Gertrudis, Tex. 

Pensacola, Fla. 

Pensacola, Fla. 

Mean of 2 ships and Dry 
Tortugas, Fla. : 

Corpus Christi, Tex. 

Houma, La. 

Houston, Tex. 

Tarpon Springs, Fla. 

Hatteras, N. C. 

Nantucket, Mass. 

Nr. Dry Tortugas, Fla. (3) 

Miami, Fla. 

Wilmington, N. C. 

Merritt Island, Fla. 

Houma, La. 

Miami, Fla. 

Perdido Beach, Ala. 

Key West, Fla. 

West Palm Beach, Fla. 

Pt. O'Conner, Tex. 

Long Key, Fla. 


Bs. a 
2s 1s 2 
? 13 
21 1 
27 ll 23 
18 
10 17 
14 45 
: l 2 23 
12 16 
ll 43 a 
22 10 7 
Me 62 ll 45 
27 8 6 
me 29 «(ll 11 
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*Sept.30, 
Aug. 13, 
Aug. 4, 
Aug. 23, 
Sept. 4, 
Sept. 5, 
Sept.16, 
June 16, 
Sept. 2. 

"Sept. 4, 
Nov. 4, 
July 3l, 
Sept.18, 
Sept.2l, 
Aug. 7, 
Aug. ll, 
Sept.23, 
Oct. 7, 
Aug. 30, 
July 27, 
Sept.14, 

*Sept.14, 
Oct. 18, 

*Oct. 19, 
Aug. 27, 
Sept.15, 
Sept.17, 

*Sept.19, 
Oct. 15, 
Sept.2l, 

*Sept.22, 
Oct. 5, 
Aug. 24, 
Aug. 26, 
Oct. 3, 
Aug. 30, 
Aug. 30, 
*Aug. 31, 
Sept.ll, 
Cct. 15, 
Aug. 12, 
Sept.19, 
Sept.23, 


Table I - Continued 


December, 1958 


Be 

DATE Po | Po Po Po Po | 
in. in. mb. mb. mb. a >) i 
1929 28.80 29.96 1.16 975 1015 40 58 65 i 
1932 27.83 30.11 2.28 942 1020 78 12 
1933 28.80 29.96 1.16 975 1015 40 24 2 iz 
1933 28.63 29.48 970 998 28 54 
1933 27.98 29.98 2.00 948 1015 67 
1933 28.02 30.24 2.22 949 1024 75 30 20 ia 
1933 28.25 29.82 1.57 957 1010 53 42 
1934 28.52 29.94 1.42 966 1014 48 37 ae 

1935 26.35 29.92 3.57 892 1013 121 6 
1935 28.71 29.89 1.18 972 1012 40 51 ‘ae 
1935 28.73 - 973. - - 
1936 28.46 30.00 1.54 964 1016 52 19 a 
1936 28.53 29.42 89 966 996 30 34 
1938 27.86 29.52 1.66 943 1000 57 50 ce 
1940 28.76 29.75 99 974 1008 34 11 zt 
1940 28.78 30.02 1.24 975 1017 42 26 27 
1941 28.31 29.66 1.35 959 1004 45 21 oF 
1941 28.98 30.19 1.21 981 1022 41 i8 ae 
1942 28.07 29.64 1.57 851 1004 53 18 e 
1943 28.78 30.02 1.24 975 1017 42 17 ls i 
1944 27.88 30.66 2.78 944 1038 94 49 17 
1944 28.31 29.39 1.08 959 995 36 26 2a im 
1944 28.02 29.80 1.78 949 1009 60 27 ia 
1944 28.42 29.67 1.25 962 1004 42 34 3 
1945 28.57 30.13 1.56 968 1020 52 18 ae 
1945 28.09 30.00 1.91 951 1016 65 12 24 
1947 27.76 29.83 2.08 940 1010 70 1g 
1947, 28.53 29.70 81.17 966 1006 40 28 33 
1947 28.59 29.65 1.06 968 1004 36 13 ‘= 
1948 27.62 29.61 1.99 935 1003 68 7 as 
1948 28.41 29.83 1.42 962 1010 43 16 a 
1948 28.85 29.77 -92 977 1008 31 27 2 
1949 28.86 30.20 1.34 977 1023 46 24 at 
1949 28.16 30.12 1.96 954 1020 66 22 23 
1949 28.45 29.95 1.50 963 1014 51 is 20 
1950 28.92 29.71 979 1006 27 21 
1954 28.35 - 960 - 
1954 28.38 « - 961 - 45 22 
1954 27.97 29.77 1.80 947 1008 61 37 32 
1954 27.66 29.55 1.89 937 1001 64 18 so 
1955 28.40 29.77 1.37 962 1008 46 45 — 
1955 28.51 29.87 1.36 966 1012 46 50 the 
‘ 1956 28.76 29.83 1.07 974 1010 36 29 so 
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28.80 975 Panama City, Flea. 
27.83 942 East Columbia, Tex. 
28.98 981 Brownsville, Tex. 
28.66 971 Cape Henry, Va. 
27.98 948 Jupiter, Fla. 

28.07 951 Brownsville, Tex. 
28.25 957 Hatteras, N.C. 

966 Jeanerette, La. 
26.35 892 Long Key, Fla. 

29.18 988 Columbia, S. C. 

238.73 973 Miami, Fla. 
28.46 S64 Fe.Walton, Fla. 
232.52 966 Mean of 2 ships off Hatteras, N.C. 
28.04 950 Hartford, Conn. 
28.87 9373 Port Arthur, Tex. 
23.78 Y¥75 Savannah WSO, Ga. 
28.66 971 Houston WBO, Tex. 
29.00 982 Carrabelle, Fla. 

28.10 952 Seadrift, Tex. 
28.78 975 Ellington Fld., Tex. 
27.97 947 Hatteras, N.C. 
28.31 956 Pt. Judith, R.I. 
28.02 949 Dry Tortugss, Fla. 
28.42 962 Sarasota, Fla. 
28.57 968 Pslacios, Tex. 
28.09 951 Homestead, Fla. (FECRR) 
27.97 947 Hillsboro, Fla. 
28.57 968 New Orleans, WBO, La. 
28.75 974 Savannah, Ga., WBAS 
28.45 963 Key vest wWBO, Fla. 
28.47 964 Clewiston, Fla. 
28.92 979 Miami, Fla. 
28.86 977 Diansond Shoals L/S, N.C. 
28.17 954 W. Palm Beach WBAS, Fla. 
28.88 973 5 mi. west Freeport, Tex. 
28.92 3979 Ft. Morgan, Ala. 
28.35 960 Aircraft Recon. 
25.42 962 Suffolk County AFB, RI- 
27.97 947 40°N = 71°W Aircraft Recon. 
27.70 3938 Little River, N.C. 
28.40 962 Ft. Macon, N. C. 
28.63 970 Cherry Point, N. C. 
28.78 975 Destin, Fla. 
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Note: All values are estimated except py. 


Po 


Pa 


Station - at 


HY 7 December, 1958 


Footnotes for Table I 


central pressure 

asymptotic pressure 

pressure difference between “outside” of storm and center 
radius to region of maximum wind speed 


(a) computed from pressure profile 
(b) observed from wind speed record 


forward speed of the storm 
time interval over which storm movement is averaged to obtain c 


lowest pressure detected by a barometer 
which pa was observed 


minimum distance from station to track of storm center 


same hurricane as previous line, entering or passing coast at a 
second point 


Lowest in region where pressure profile parameter were computed 
29.31 inches observed at Mobile, Ala. later 


Parameters obtained by interpolating is time between ship off 
western end of Cuba and Miami, Fla. and apply to vicinity of Dry 
Tortugas. 
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Fig. 2. Reconstructed surface pressure, wind speed and deflection angle 
for 1200 EST in the hurricane of September 21, 1938, 
After Myers and Jordan. 
(A) Surface Pressure, (inches of Mercury) 
(B) Wind Speed (M.P.H.). X is location of pressure center 
(C) Deflection Angle (degrees) 
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Eq. (4) it is assumed that the winds will blow parallel to the isobars. A more 
sophisticated theory, as well as observations show that the wind must cross 
the isobars toward lower pressure. The angle of inflow, that is the angle with 
which the wind crosses the isobars, must be asymmetric with respect to the 
center of a moving storm. The best determination of the angle of inflow for 
New England hurricanes is shown as Fig. 2c. (8) 


Tides and Storm Surges 


If the effects of wind waves and swell are eliminated, the resulting ele- 
vation of the sea surface is usually controlled by the gravitational attraction 
of the sun and moon and the topography of the region of observation. Minor 
changes of a seasonal nature are also produced by the seasonal changes in 
atmospheric pressure, prevailing winds, and temperature of the water. The 
effects of these gravitational and climatic factors can be computed many 
months in advance and in this paper are identified as the normal tide. 

In the neighborhood of an intense storm, the actual level of the sea will 
differ significantly from the normal tide, and the actual level of the sea may 
be referred to as a storm tide. The difference between the storm tide and the 
normal tide is called the storm surge. 

Economically important storm surges may be produced by extratropical 
storms as well as hurricanes. Extratropical storms have been responsible 
for the highest storm tides north of Cape Cod and at a few more southerly lo- 
cations. However, the surge produced by many hurricanes is much higher 
than those produced by any extratropical storm. Fig. 3 shows the maximum 
recorded tide at a representative number of places along the Atlantic and Gulf 
Coasts of the United States.(9) The remainder of this paper will be devoted 
to tropical cyclones and the storm tides produced by them. 

The peak water level reached in a hurricane will always be somewhat 
greater than the storm tide level because of the effect of waves with periods 
of only a few seconds. It is desirable to separate these two phenomena be- 
cause their effects are quite different. If the peak tide level is two feet above 
mean sea level, it will penetrate inland to the contour which is approximately 
two feet above mean sea level, even if this distance is several miles, and it 
will not rise appreciably above two feet regardless of the slope of the beach. 

A two foot wave, on the other hand, may run up to a sloping beach to a 
height of seven or eight feet above the still water line, but under most con- 
ditions the effects of such a wave cannot penetrate inland more than a wave 
length from the still water line, and most of the damage from such waves will 
be confined to a much narrower trip. The waves which accompany the land- 
fall of a hurricane will generally be in the range of 100 to 1000 feet in length. 
The remainder of this paper will be restricted to a discussion of the storm 
tide. 

In most cases the principal cause of the storm surge is the frictional drag 
between the wind and water. Changes in atmospheric pressure, although usu- 
ally of secondary importance, may occasionally be more important than the 
wind effect.(10) Waves breaking against the shore or over barrier reefs may 
also contribute to the piling up of water near the shore. 

Fig. 4 shows the storm surge produced by Hurricane CAROL at 19 Coast 
and Geodetic Survey tide stations. The dots on the curves for New London, 
Newport and Woods Hole indicate approximate surge heights after the tide 
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Fig. 3. High tide data for stations having a record for 20-years or longer. 


gages became inoperative due to the storm. The dashed line crossing the 
storm surge curves indicates the time at which the storm was nearest each 
tide station. The storm track is shown to the left of the storm surge curves. 
The letters on the map and beside the curves identify the storm surge curves 
with the geographical location of the tide gage. The maximum surge height 
and mean tide range are shown on the left. 

With the exception of Willets Point, near the western end of Long Island 
Sound, the maximum surge height at all stations between Hampton Roads and 
Cape Cod occurred between the passage of the storm and two hours later. At 
Willets Point, the maximum surge occurred more than four hours after the 
passage of the storm center. Thus, the water continued to rise at many places 
to the left of the storm track after the wind had acquired an offshore com- 
ponent, suggesting that the surge arrives as a wave which had been generated 
elsewhere. This is particularly true at Willets Point. The notion of a wave 
receives further support from the damped oscillations readily seen in the 
records for Atlaatic City, Sandy Hook, and the Battery, and suggested in the 
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records of several other stations. These oscillations, calied resurgences by 
Redfield, (4) are characteristic of the hurricane surge in this region. 

Fig. 5 shows the storm surge as a function of time at the recording tide 4 
station nearest the point of maximum storm surge height for several hurri- 
canes. A trend toward rising or falling sea level may become evident 12 to D's 
24-hours before the arrival of the storm. In either, a rapid rise begins when -_ 
the storm is about 100 to 200 miles distant, that is when the wind speed reach- 
es about 30 miles per hour with the approach of the storm. This period of & 
rapid rise lasts from two to six hours and is generally followed by a slightly q 


more rapid drop to something below the normal water level in areas with good E 
drainage. In marshland or other areas with poor drainage, the water may re- ‘ 
main well above normal for several days. Fig. 6 shows the contours of the a | 
storm surge along the coast for four hurricanes. Notice the tendency for the - 


peak value to be higher on the right or to extend farther from the center on 
the right hand side of the storm track. Notice also that the peak storm tide 


may be as much as 10 feet higher than the storm tide 60 miles to either side 
of the peak. 
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Fig. 7 based on data furnished by the Corps of Engineers, Coast and Geo- 
detic Survey, Woods Hole Oceanographic Institution and others, shows the 
storm surge heights produced by Hurricane CAROL at a great many locations 
along the southern New England and Long Island coasts. Here the effects of 
topography and local wind conditions are shown by the higher surge heights 
near the heads of converging bays and at exposed points, and the decreased 
surge heights near the widest section of Long Island Sound. The varied effects 
of topography and local wind conditions make it difficult to construct graphs 
such as Fig. 6, or to develop any reiation between storm surge height and 
meteorological parameters. It is clear that one should avoid using surge 
heights obtained near the heads of estuaries in developing such relationships, 
but it is not possible to avoid all local effects. 


Maximum Storm Tide as a Function of Storm Intensity 


The effect of wind in piling up water is nearly proportional to the wind 
stress. The wind stress is generally taken as being proportional to the square 
of the wind speed, and according to Eqs. (2) and (4), this is approximately pro- 
portional to the pressure deficiency in a hurricane. The direct hydrostatic ef- 
fect of the decreased pressure is also proportional to the decrease. Thus, we 
are encouraged to look for a correlation between (Py - Po), or Po alone, and 
the maximum storm surge, that is the maximum effect of the storm on the 
level of the sea. The results of such an effort are shown in Fig. 8. 

The highest observed tides along the Gulf Coast of the United States during 
30 hurricanes is indicated by dots as a function of the central pressure of the 
storm as it crossed the coast. An attempt has not been made to eliminate the 
normal tide, as the normal tide range is less than two feet in most of this 
region, and sufficient data for the elimination of the normal tide were not 
available for all of the reports. The regression equation based on these data 
gives a correlation coefficient of 0.68. The extreme differences between the 
storm and the normal tide for 10 Atlantic Coast hurricanes are indicated by 
x’s. The tide range along the Atlantic Coast is much greater than in the Gulf 
of Mexico, and a correction for the normal tide is necessary in order to ob- 
tain any useful degree of correlation. Largely because the variability due to 
the normal tide is eliminated, the correlation coefficient for Atlantic data has 
been increased to 0.72 even though these data were not considered in the deri- 
vation. Consideration of this single parameter is thus sufficient to account 
for half of the observed variability in maximum storm surge heights. (12) No 
improvement was obtained by considering P, and R as defined in Eqs. (3) and 
(4), or by using the storm tide potential as defined by Reid. (6) 


Effects of the Forward Motion of the Storm 


Proudman(13) using a very elementary theory, has shown that the response 


of the water in a rectangular canal to a moving atmospheric disturbance could 
be expressed in the form: 
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Fig. 8. Storm surge height, as a function of the central pressure 


in the storm. 
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height of water level disturbance 


Atmospheric pressure disturbance expressed as the height of a 
water barometer 


speed of atmospheric disturbance 
depth of canal 
= acceleration of gravity 


Later writers(14) have shown that these results are modified somewhat, when 
the effects of the earth’s rotation and the two dimensional nature of the dis- 
turbance are considered. The essential feature of the development which 
carries over to more complex situations is, that for any basin there is a criti- 
cal velocity of motion of the atmospheric disturbance which will produce maxi- 
mum response in the water. An examination of the residuals obtained from 
Fig. & supports this concept, but there are not sufficient data from any type of 
coastline to establish the proper form of this relationship from empirical 
evidence. 

The speed of the storm can become important in quite another way, for the 
speed of the storm, together with the radius of maximum winds determines 
the duration time and the length of the fetch over which the storm winds can 
act to produce storm tides. 


Hurricane Frequency 


From a practical point of view, it is necessary to consider the frequency of 


storms as a function of intensity, size, and region of occurrence. Any real 
trends in such frequencies would also be highly significant in an economic 
sense. Unfortunately, the amount of data available for making such studies is 
rather limited and all data are not of the same quality. The nature and sense 
of the available data can be indicated by three figures adapted from recent 
Weather Bureau reports to Congress, (16) 

Fig. 9 shows the frequencies of tropical storms and hurricanes for the 
period 1887-1956. A tendency for years of above normal hurricane frequency 
to be followed by a period of below normal hurricane frequency can be easily 
noted. There is no evidence of any well marked trend or regular cycle in 
hurricane frequency. It may be possible to obtain a useful indication of the 
number of storms to be expected during the next two or three years by con- 
sidering the number which have occurred during the last 10 to 20 years, but 
the best estimate of hurricane probabilities more than a very few years in the 
future would be obtained from considering the frequencies for the entire peri- 
od of reliable records. 

Fig. 10 shows the number of times destruction has been caused by tropical 
storms on the United States mainland during the present century. 

Fig. 11 shows the accumulated freauencies of tropical storms, whose 
central pressures were 29.00 inches of mercury or less, for several sections 
of the U. S. Coast. Cape Hatteras is taken as the dividing line between the 
North Atlantic and South Atlantic sections of the coast for this figure. The 
data for the Florida Keys are clearly different from those of the other regions. 
Neither the size of the sample available, nor the present state of theoretical 
knowledge concerning hurricanes is sufficient to determine whether the 
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Fig. 9. Chronology of storm tracks. The occurrence of hurricanes and 
tropical storms in the different coastal regions shown at 
left, and the total occurrence, 
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Fig. 10. Number of times destruction was caused by tropical storms, 
1901-1955. 


differences in the other five curves are due to sampling error or to some 
fundamental aspect of hurricane behavior. 


Hurricane Observations and Forecasts 


If a hurricane approaches land in the neighborhood of a powerful weather 
radar installation, such as that now operating at Cape Hatteras, N. C. and on 
order for about 15 other coastal sites, an experienced radar hurricane ob- 
server can locate the center of the storm within 5 to 10 miles while the storm 
center is within 150 miles of the station. The accuracy decreases as the 
distance from the radar station is increased. 

When the storm is out of range of land based radar, the position of the 
storm center must be determined from the available ship and aircraft reports. 
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Occasionally the eye of the storm will pass directly over a ship. Somewhat 
more frequently, but not on every day of the hurricane’s existence, an airplane 
will penetrate the eye of the hurricane. In these cases the storm center can 
often be located within an accuracy of about 10 miles. In the more usual case, 
the storm center must be located in the center of an isobaric system de- 
termined from fewer than a dozen ship reports, or from radar observations 
made from aircraft on the edge of the storm. The error in locating the storm 
center in this way will frequently be as much as 20 or 30 miles and occasion- 
ally as much as 60 miles. Errors in communication via radio and teletype- 
writer compound the uncertainties resulting from the currently unavoidable 
shortcomings in our observational procedures. 

A few months after the storm, when mailed in reports, free of communi- 
cation errors, as well as much additional information not available to the fore- 
caster during the storm are obtained, it is possible to repeat the analysis and 
determine the position of the storm track with an accuracy impossible during 
the storm. 

A simple extrapolation of the storm path accounts for considerably more 
than half of the obtainable skill in hurricane forecasting. Even the more 
sophisticated forecast methods, based on a knowledge of flow patterns at all 
levels, give only a velocity and acceleration of the storm movement, and per- 
petuates any existing error in the storm position. 

When all of these factors are considered, it is evident that storm 
movements cannot be forecast with timetable accuracy. A study of forecast 
accuracy for the period 1952 to 1955 indicated an average error of 69 nautical 
miles in the 12-hour forecast of storm pesition.(15) This was made up of an 
average error of 83 miles during the period 1952-1953 and 63 miles in 1954- 
1955. A gradual improvement in forecast accuracy is to be expected as the 
result of improved instrument and forecast techniques. 

A hurricane watch is announced for specific areas whenever it is recognized 
that a specific hurricane threat exists for any region of the coast. The hurri- 
cane watch is announced before it is possible to delineate the area in which 
the storm is likely to go inland and when the hurricane effects are believed to 
be at least 24 hours away from the coast. A hurricane warning is issued for 
a specific area whenever there exists a reasonable probability that the storm 
will go inland in that region within the next 24 hours. When the average accu- 
racy of forecasts and the length of coastline seriously affected by a hurricane 
are considered, it is evident that the hurricane warnings must cover an area 


several times as large as that which will actually suffer severe damage from 
the storm. 


SUMMARY 


A full realization of the hurricane’s potential for danger and of the difficulty 
of predicting its future behavior exactly, should materially reduce the loss of 
life and property due to these storms. 

More than 390(1) persons lost their lives in Hurricane AUDREY, largely 
because they had lived through other hurricanes with comparatively little 
damage and did not believe that this storm could be as severe as predicted. 

Hurricanes vary in intensity. Tropical storms which barely qualify as 
hurricanes may elevate the sea level no more than four feet, but severe hurri- 
canes may elevate the sea level by 15 feet or more above normal tide for the 
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time of the hurricane landfall. Intense extratropical storms may produce 
rises of five to six feet above normal. 

The sea level may be raised by one or two feet for a thousand miles of 
coastline by even a moderately intense hurricane, but the really severe storm 
surges of more than 5 feet above the normal tide level rarely extend over 
more than 100 miles of coastline. 

The duration of the rapid rise in sea level is rarely greater than 6 hours, . 
and in areas with good drainage the fall is even more rapid. In areas of poor 
drainage, moderately high water may persist for days after the storm. 

Knowledge of the intensity of a hurricane before it reaches the coast is 
meager. The average error in a 12-hour forecast is approximately 60 nauti- 
cal miles, and if warnings are to be effective, an area several times as large 
as that actually affected must be alerted. 

Disasters such as that which occurred with Hurricane AUDREY can be 
eliminated by the evacuation of the coastal plain in the threatened region to an 
elevation of 10 to 20 feet above the normal high water level. In marshland, 
where high ground may be as much as 40 to 100 miles away by highway, many 
mass evacuations may be required to avoid one catastrophe. With improved 
instrumentation and forecasting techniques, the number of evacuations, later 
proved to be unnecessary, will be reduced, but no technique now available will 
permit the prediction of hurricane motion with timetable accuracy. 
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MECHANICS OF SEDIMENT-RIPPLE FORMATION#@ 


Corrections by J. Bogardi 


J. BOGARDI.!—Corrections to Discussion by J. Bogardi: 


Page 1832-3, 14th line “. . . difficulties to Mr. Liu—of the multitude. . .” a 
should be “difficulties to Mr. Liu to determine these correlations—of the a 


multitude.” 
Page 1832-4, Eq. (8) “T, = U*2 = ” should be “T,, = pU*= ”. 


Corrections to Closure by H. K. Liu: 


Page 1832-29, 16th line “. . . This factor is proportional to gd/U,?2 ieee # 
should be “The third parameter f is proportional to d/DS = gd/U,2, the ¥ 
channel stability factor, which is the inverse... .”. 

Page 1832-29, delete (a) the words “in the form” before Eq.(35), (b) Eq. 
(35), and (c) the line following Eq. (35), and substitute the following: 

“If one uses the channel stability factor d/DS = gd/U,2 the configuration of 
the bed is quantitatively defined by the parameter 8 mentioned above unique- 
ly. The relationship of B and gd/U,2 according to Bogardi is 


= (35) 


Y 7 
From the data of Fig. A, Bogardi assumes that N = 0,882. 


Bogardi further states: Obviously there is a parameter uw proportional to 
U,d/v , which defines the bed configuration uniquely too and 


(36) 


Finally Bogardi pointed out that a third parameter € proportional to U, 
naturally defines the bed configuration uniquely too: 


(37) 
B, a and € are obviously interrelated, and their value depends upon the 
temperature, as well as on the specific gravity of the sediment. Assuming a 
temperature of 20° C., and a specific gravity of 2,65, the parameter values 
are obtainable directly from Eqs. (35), (36), and (37). 


Fig. 19 shows after Bogardi the U, = eas equations. 


. Proc. Paper 1197, April, 1957, by H. K. Liu. 
1. Asst. Prof. Technical Univ. of Budapest, Hungary. 
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According to Fig. 19 ripples are only possible if 1,777 < € < 4,06 and 
d< 2mm. Bogardi constructed similar figures for Eqs. (35) and (36).” 
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NORTHEASTERN FLOODS OF 1955: RAINFALL AND RUNOFF®@ 


HOWARD M. TURNER, ! M. ASCE.—This is a very interesting engineering 
description of the floods of August and October, 1955. The data selected and 
the presentation of it give an excellent summary of the magnitude and extent 
of this storm, which, as the author points out, exceeded previous records of 
intensity of rainfall and flood flows and covered a tremendously wide area. 
The writer has made some further comparisons of the magnitude of this flood ie 
with other large floods in New England. eae 

Fig. 1 shows area-depth curves for total rainfall in various New England 
storms showing how the 1955 storm exceeded previous records even of the 
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very large four-day hurricane storm from Septen:ber 17-21, 1938. Most of 
the 1955 storm came the two days of August 18 and 19. 

The Storm Rainfall report of the Miami Conservancy District, Revised, 
1936,(1) gives area-depth curves of great storms in the northeastern part of 
the country from 1869 to 1933. On Fig. 2 the 1938 and 1955 hurricane storms 
are plotted on the two-day area-depth curves for the Northeast taken from 
that report showing how much the 1955 storm exceeded any of these record 
storms. 

It is interesting to compare past analyses of storm intensities with the 
later records and with the 1955 storm and to note how the results have 
changed by records which have been added through the years. As an example, 
Fig. 3 shows a curve of high rainfall frequency taken from “The Frequency 
of High Rates of Rainfall” by Allen Ilazen, published in 1921.(2) The writer 
has plotted on this David L. Yarnell’s values, published in 1935(3) for 10- 
year and 100-year maxima estimated for Massachusetts and also the maxi- 
mum actual rainfalls for different periods of time recorded during the 1955 
storm. The results of the latter would require an entire recasting of these 
frequencies derived thirty-six years ago. 

In a paper of “Flood Flow Formulas based on Drainage Basin Character- 
istics” by H. B. Kinnison and B. R. Colby, (4) a frequency curve of storm rain- 
fall was included. This is shown on Fig. 4. On this are plotted the maximum 
records for the 1955 storm showing the extreme magnitude of this storm com- 
pared with that curve. 

It is naturally expected that a longer period of records would have the re- 
sult of increasing the maximum results obtained, but it is interesting to note 
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that the present period does not seem to show any orderly increase as a long- 
er time of records alone would cause. Fig. 5 shows records of large three- 
day storm rainfall since 1855. The two largest storms have come in the last 
19 years, exceeding by 25 per cent anything of record before. 

Considering flood flows, records of the gage heights on the lower 
Connecticut River at Hartford are available for the last 114 years. There are 
also a few high-gage records going back for 318 years. The greatest flood 
was in 1936; the next, in 1938; and now the third, in 1955. As far as can be 
determined from storm and flood histories, nothing of these recent large 
magnitudes occurred in over 300 years. Fig. 6 shows this large flood record 
of the Connecticut River at Hartford back to 1683. From all that can be 
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determined, the three largest floods have occurred in the last twenty years. 
The great size of the 1955 storm is shown by the large flood it caused on this 
river though it was largely confined to the lower third of the drainage area. 

The author has shown the relation between the peak flow and drainage 
areas on his Table 5 and Fig. 6. The writer has made up Fig. 7 showing the 
the same data for New England, including the large floods previous to that of 
August, 1955. A curve of the 40% Myers rating is shown as an enveloping 
curve including practically all records up to that time. The 1955 flood re- 
quires a 70% rating curve to include the high records of the 1955 flood for the 
drainage areas of 300 to 700 square miles. 
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These very high rates of flow on these fairly small drainage areas shown 
in the Myers rating and in the rainfall frequency studies are due to the heavy 
rainfall in a short time. The author shows on Fig. 3 the rainfall curve at 
Hartford with the total of 14.4 inches in a total time of about 42 hours. Actu- 
ally 97 per cent of it fell in 28 hours; and about 83 per cent of it, in 24 hours. 


Runoff 


The author gives some figures of the percentage of runoff to rainfall, 
which, except for a few cases, appear to be somewhat lower than might be ex- 
pected from a flood of this intensity. Some of the very high ones that are 
given may very likely have been affected by the failure of dams. 

The author has noted the rapid rise and fall of the hydrographs of this 
flood. In the case of the large rivers the reason for this is quite evident. 

For example, on the Connecticut River the heavy storm flows were confined 
to the lower third of the drainage area. Channel storage was thus confined to 
only the lower portion of the river valley, and the recession curve of the 
hydrograph representing the drainage of the channel storage would be shorter 
in duration as the amount of storage to drain away would not be so great. The 
same thing shows in the recession limb of the hydrograph of the Delaware 
River at Port Jervis, where the heavy storm rainfall was concentrated into 
about 50 per cent of the drainage area. 

The hydrographs of the smaller rivers remain to be analyzed. The writer 
believes, however, that the 6- or 12-hour unit hydrographs for these rivers 
will prove to be close to the corresponding unit hydrographs for the 1938 
flood except that the peaks may be somewhat higher. 


Frequency 


The question of frequency is one of the most difficult problems in flood 
studies, but some frequency rating of the size of a given flood seems to be re- 
quired. The author gives a tabulation and a frequency curve for the North- 
east, according to the United States Geological Survey’s method of rating the 
recurrence interval of any given flood as its ratio to the mean annual flood, 
but the curve gives the recurrence interval up to 100 years only, not long 
enough for many requirements. This curve shows very different results from 
that developed for Connecticut by B. L. Bigwood and M. P. Thomas. (5) 

The section on flood frequency of the 1942 report of the Committee on 
Floods of the Boston Society of Civil Engineers, published in 1942, (6) analyzed 
the frequency of floods on the Connecticut River at Thompsonville, where 
there were then over 93 years of record. Various methods were used. The 
results, using data up to 1935, i.e., prior to the 1936 and 1938 floods, gave a 
frequency of the 1936 flood from 10,000 to 16,7000 years, depending on the 
method used. The same computations made with three years’ more data in- 
cluding these two large floods gave frequencies from 300 to 400 years. 

Apparently the same difficulty is still found by the United States Corps of 
Engineers as shown in its recent report on the 1955 flood,(7) where for four 
rivers in southern New England the flow for a given frequency used in 1956 
was much higher than shown by the curves tentatively adopted prior to the 
1955 flood. See Table 1. 

The paper by Kinnison and Colby, published in 1945, (4) which has been used 
a great deal in Massachusetts, rates floods as “minor,” “major,” and “rare” 
with stated frequencies of 15, 100, and 1,000 years respectively. On the rivers 
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on which there are records concerning the 1955 flood, the Kinnison-Colby 
figures would be as shown on Table 1. These generally rate the higher 1955 
flood flows as equal to or exceeding a “rare” flood. 

In that paper the flood frequency ratings were based upon precipitation and 
depth of runoff. The use of precipitation in determining flood frequency has 
its difficulties because of the varying percentage that occurs in the form of 
runoff in the flood, but the rainfall data are more inclusive. The writer be- 
lieves that the flood runoff figures would have some advantage as data for 
flood frequency computations instead of the peak floods. The runoff data are 
often not so readily available and require more work to use, but the results 
are on a broader basis, eliminating the varying characteristics of each 
stream. With a runoff frequency established, the peak flood of each location 
could be determined by unit hydrograph methods. On a few of the streams 
for which Kinnson-Colby figures are available for the high 1955 flood areas, 
the runoff compares to their “rare” or 1,000-year flood. 

All this points out the difficulty of the problem. Each of these successive 
floods gives more information with which to work; but as the author points 
out in his paper, figures that still show such wide variation are not very 
“satisfactory.” Yet, it is a problem which is important. In over 950 years 
of records of the flow of the river Danube, the greatest flood was in 1501; the 
next in 1787 was 85 per cent of it; and the third one in 1899, 75 per cent of 
it.) There is no question but what these extraordinary events occur at wide 
intervals, and we have to know something of the degree, time, and magnitude 
in flood works. 

Fortunately, frequency is not required so much in the design for spillway 
capacities of dams except indirectly; but for other designs, such as flood- 
control reservoirs, flood channel improvements, and other flood-control 
work, some method is needed to fit the design flood to the value of the pro- 
posed construction. The 1955 flood has certainly given perhaps an unscien- 
tific but in some ways very convenient yardstick. It has been taken as the 
design flood on some recent river-improvement designs. As far as the public 
is concerned, a large storm like that of 1955 is, of course, a very good flood 
measure, It is much easier to describe flood-control measures for a flood 
intensity that has already been experienced than it is to talk about the some- 
what mythical frequency of a 100- or 1,000-year flood, which may come in the 
next year. 

It is interesting to look back and see the tremendous development of the 
science of flood analysis and prediction since the 1927 flood. There are still 


problems, but the progress made justifies the belief that the future will solve 
them. 
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EXPERIMENTS ON SELF-AERATED FLOW IN OPEN CHANNELS 
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SYNOPSIS 


Measurements of distribution of air concentration in self-aerated flows 
are presented in this paper. The experiments were made in a rough channel 
of sand-grain type surface at various slopes and discharges, and the data 
were used as a basis for study of the mechanism of entrainment of air and to 
relate the air content and distribution to the characteristics of the flow. The 
analysis of the data shows that the air distribution can be adequately de- 
scribed by relationships based upon a simplified concept of turbulent transport 
and thus are functions of the flow characteristics. The maximum depth and 


the mean velocity are both shown to increase above those of a corresponding 
nonaerated flow. 


INTRODUCTION 


A characteristic of high-velocity, open-channel flow is the phenomenon of 
self-aeration in which atmospheric air is insufflated into and mixed with the 
flow to create the appearance of “white water” with its violently agitated and 
ill-defined free surface. This condition is frequently observed in flows down 
steep chutes and spillways. 

It has been reasonably well established by several observers(1,2,3) that 
for flows over a spillway from a quiescent reservoir, incipient aeration does 
not occur on the slope until a point or region is reached at which the boundary 
layer thickness is equal to the depth. The characteristic roughening of the 
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water surface immediately upstream of the appearance of “white water” can 
be readily observed in such occurrences. With increasing initial water 
depths, the points where the surface roughening and the “white water” occur 
move downstream along with the intersection of the boundary layer with the 
water surface. These considerations suggest that the aeration phenomena are 
related to the conditions of turbulence in the flow. Examination of such flows 
by means of high-speed photography indicates that the upper boundary of the 
flow is rather ill-defined. It consists of a zone which appears white because 
of entrained air, while above this zone a spray of water droplets occurs mov- 
ing more or less parallel to the flow and below which is a region of discrete 
air bubbles suspended in the fluid. Other investigators have described air- 
entrained flow as being made up of air bubbles entrained in water in the lower 
levels and water particles in the air in the upper regions. This concept is 
further emphasized by measurements of the air concentration distri- 
bution(4,5,6) which show that the air concentration increases continuously 
from the bed, exhibiting a smooth transition from a finite value near the bed 
to a maximum value of 100 per cent at the free surface which value is ap- 
proached asymptotically. 

This paper presents the results of systematic experiments on distribution 
of air in self-aerated flows in a channel set at various slopes and operated at 
various discharges. A provisional analysis is made of the data with respect 
to air distribution in a vertical section transverse to the direction of flow 
with a view of arriving at empirical relationships between concentration 
parameters and the flow characteristics. 


Experiments and Results 


Apparatus 


The experiments were undertaken in a channel 50 ft. long, 1.5 ft wide, and 
1 ft deep, which could be adjusted to any slope from horizontal to nearly verti- 
cal. The width was chosen so that for ordinary flows sidewall effects due to 
growth of the sidewall boundary layer and air entrainment along the sides 
could be avoided in the central portion of the cross section. That the width 
was adequate to make the flow two-dimensional in the center portion was 
demonstrated by earlier measurements of air concentration and velocity in a 
transverse section at the end of the flume.(4) The flow was obtained by gravi- 
ty from the laboratory main supply flume through two feed lines and controlled 
by hydraulically operated valves. The water reached the inlet of the flume 
through hollow support members and through a rectangular conduit on the 
underside of the channel. Fig. 1 shows the channel and appurtenances. 

At the inlet the flow is guided through two 90-degree vaned bends and a 
contraction. The inlet gate itself is cantilevered from its mounting at the top 
of the flume so that there are no slots or guides to disturb the flow. Further 
guidance is provided by a lip attached to the bottom of the gate and extending 
inward to further contract the flow; therefore, the gate opening corresponds 
to the depth of flow at the upstream end of the channel. The maximum gate 
opening is 0.5 ft. Pitot measurements of the flow just downstream of the inlet 
with an 0.25-ft gate opening showed no deviation greater than 2 per cent from 
the mean velocity except in the bottom corners where slightly smaller veloci- 
ties were encountered. 

In the experiments here described, the air-entrainment process was in- 
tensified by the installation of artificial roughness on the channel bed that was 
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considerably rougher than the painted steel sidewalls of the channel itself. 
The material used for the roughness was a commercial non-slip fabric coat- 
ed with granular particles. The particles, whose mean size was 0.028 in., 
had a mean spacing of 0.039 in. and were embedded in a mastic, as shown in 
Fig. 2. ; 

The discharge rate was measured by two calibrated Pitot cylinders mount- , ioe 
ed in the supply lines and leading to differential gages mounted on the control 5 
panel. Each discharge meter was calibrated separately for flows from 1/2 
to 10 cfs with an accuracy of the order of 1 per cent. 

A control panel contained the actuating components for the hydraulic 
system and indicators for the measurement of slope, gate opening, and dis- 
charge which permitted adjustment of these variables to predetermined 
values. 

The electrical air-concentration measuring instrument has been fully de- 
scribed in an earlier report. (7) The method consists basically of a measure- 
ment of the difference between the conductivity of a mixture of air and water 
and the conductivity of the ambient water alone. A strut supporting a pair of 
probes is combined with the electrical circuit so that the air-concentration 
measurements may be made in a small region of the flow. The probes were mer 
1/4 in. wide and 1/4 in. apart. The instrument is direct reading, and the eh 
relation between conductivity and air concentration can be determined ana- a 
lytically. The strut holding the probes is arranged so that it is possible to : 
traverse the flow cross section both vertically and laterally. The instrument Se, 
is limited to measurements not closer than 0.02 ft from the walls because of % an 
current deflections by the walls at lesser distances. Since the measurements 4 
of the air concentration depend upon the resistance or specific sensitivity of 4 
the suspending water, the meter was zeroed in water continually circulated 
from the test channel through a deaerating tank, thus providing water of the aes 
same character as that in the flume. Oscilloscope measurements show a 
rapidly fluctuating current between the probes. Since the meter could not 
follow these rapid fluctuations, the value read on the meter represented an 
average value. This average compared well with results obtained by means Fk 
of a mechanical sampler which had been previously developed at the Labora- 3 * 
tory and which was not subject to the same averaging problems. A photograph me: 
of the probes and strut supports for the air-concentration meter is shown in 
Fig. 3. 

Some of the experimental runs also required taking velocity traverses in a = 
verticals through the aerated flows. The instrumentation for making these 4 x 
observations has been previously described(8) in some detail. It consists aes 
essentially of a device for timing the travel of minute salt-water cloudlets in- sale e 
jected into the air-water mixture. The operating principle of the velocity wal 
meter consists basically of marking a small element of flowing mixture and 
then recording the time interval required for this marked element to traverse 
a fixed distance. The marking is accomplished by making a diminutive portion 
of flowing water more highly conductive to electrical current by injecting a 
small amount of salt solution into the stream. The passage of this ionized - 
cloudlet is then detected by electrodes at stations fixed 3 in. apart in the flow YF 
path. The injections were made at a speed of about 20 per sec with individual i 
injections of about 0.03 cu cm of 6-per cent saline solution injected in about 
1/600 sec. The accuracy of the air-concentration meter and the velocity 
meter was checked by comparing the measured water rate entering the 
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channel with the integrated flow determined by the velocity traverses and the 
air-concentration traverses over the channel cross section. 


Experimental Data 


For the air-distribution measurements described in this report, the 
channel slope and water discharge were varied independently. Experiments 
were made for slopes of 7.5, 15, 22.5, 30, 37.5, 45, 60, and 75 degrees. For 
each slope, measurements of air concentration were made for total-water 
discharges from 2.2 to 9.6 cfs and for some slopes up to 15.0 cfs. 

The measurements were made at Station 45 (45 ft from the inlet) in a verti- 
cal centerline plane. Readings were taken at 0.01-ft intervals on a line 
normal to the channel bottom with the lowest point 0.02 ft from the bed and 
the uppermost point where the probes were completely out of the flow and the 
meter registered zero water concentration. The concentration profile was 
obtained for each slope and each water discharge after the channel had been 
adjusted for equilibrium conditions. For this purpose equilibrium flow was 
considered to be that condition of the flow of the air-water mixture for which 
the air-concentration distribution was the same at two sections 10 ft apart 
along the channel. This condition could be obtained by the adjustment of the 
initial depth of flow through control of the inlet-gate opening. Repeated 
measurements of the profiles were made at the two sections for different gate 
openings until the two air-concentration profiles were similar. Comparisons 
were made of the concentration at corresponding distances normal to the bed. 
A typical comparison of the profiles at the two stations for equilibrium flow 
is shown in Fig. 4. 

The results of the experiments in the form of measured concentrations for 
each slope and discharge have been tabulated in Appendix B for reference. A 
tabulation of computed parameters is given in Appendix C. 

The air concentration C is defined as the ratio of the volume of air per 
unit volume of air and water as measured by the concentration-meter probes. 
It is assumed that the concentration as measured represents the average 
value of the concentration in the region described by the probes of the concen- 
tration meter and is applicable to the midpoint of the probes. Essentially 
continuous curves of air concentration with respect to distance normal to the 
bed were obtained for each experiment. The concentration increases gradual- 
ly from the bed—more rapidly in the central portion and then more slowly in 
the upper region, apparently asymptotically approaching the limit of 100 per 
cent air concentration, as is shown in Figs. 4, 8, and 9. It appears that the 
complete curve is comprised of two parts which have basically different 
characteristics. Although there is no sharp demarcation between the two 
parts, the curves tend to support the description that in the lower regions of 
the stream air bubbles are suspended in water, while the upper regions con- 
sist of water droplets in air. The existence of these water droplets can readi- 
ly be detected by holding one’s hand just above the main flow. 

Since there is no definite upper boundary to the air-entrained flow in an 
open channel, questions of the definition of the linear dimensions of the flow 
arise. A number of different depth parameters related to the air concentration 
and distribution may be defined, each having particular applications in de- 


scribing the bulk-flow conditions. These depth parameters may be described 
as follows: 
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(1) A depth d is used to represent a mean depth of flow that would exist if 
all of the entrained air were removed up to the highest point where 
water is found. It therefore corresponds to the depth of a nonaerated 
flow of a given discharge with a velocity equal to that of the entrained 


flow. It may be defined as a om 

d = -—C)dy (1) 

4 


where y is the normal distance from the bed and C is the concentration 
of air as a function of y. The value of d can be approximated with suf- 
ficient accuracy by a summation of the products of the water concen- 
trations and their corresponding incremental depths. 

(2) An upper boundary of the air-entrained flow or an upper depth d, may 
be defined as the value of y where the air concentration has some arbi- 
trarily prescribed value. For some purposes, such as determining the 


mean depth, and because of the asymptotic nature of the curve, it is de 
convenient to define dy as that value of y where the air concentration C 3 

is 0.99. Other depths may be similarly defined. Such a depth has be 
previously been defined (4) as that corresponding to that value of y Sa — 


where the air concentration is C = 0.95. This depth was chosen arbi- 
trarily as a depth where the concentration of air can be reliably 
measured and which includes from 98 to 99 per cent of all of the water 4 
and consequently represents a maximum distance from the bed where a 
water will be found. 

(3) A third depth parameter may be defined on the basis of the air- 
concentration distribution as that depth dy which represents the value 
of y where the transition from the distribution in the lower regions to 
that in the upper regions occurs. Its value would then depend upon the 
characteristics of the channel as well as the air distribution. This 
depth parameter will be discussed in detail in a later section. 


The mean air concentration in the vertical can also be defined in different 
ways depending upon the particular application. The mean concentration over 
the whole range of concentrations measured is defined as 


| 
¢g (2) 
dy J y 


where C is the mean concentration of the entire flow and dy is the so-called 
upper limit where the concentration is 0.99. Another mean concentration Cy 
may be defined as the mean concentration in the region below the transition 


depth y = dy which applies to that air which is being transported by the flow, 
that is 
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Dimensional plots of d, d,,, and dy for the various discharges and for each 
slope are presented in Fig. 5 to indicate the range of values as determined 
from the concentration data. In Fig. 5(c) the mean depth d varies with dis- 
charge and slope in much the same manner as nonaerated flow. The mean 
depth increases regularly with the discharge and decreases with the slope. 
The data for the upper limit dy (Fig. 5(a)) indicate that for a given discharge 
d,, tends to decrease at first and then to increase as the slope increases; 
consequently, the maximum depth tends to be nearly the same for correspond- 
ing discharge rates on the various lower slopes. It appears that the increased 
depth due to aeration increases with both slope and discharge; however, in 
combination with the mean depth, which decreases with slope and increases 
with discharge, the upper limit tends to decrease and then increase as the 
slope increases. As a result, over a considerable range of slopes (15 to 30 
degrees), the upper limit is essentially the same and depends only on dis- 
charge. A similar situation exists in regard to dy (Fig. 5(b)) in that it also 
shows a tendency to decrease gradually and then to increase gradually as the 
slope increases. Consequently, the range of variation of dy for any given dis- 
charge is relatively little over a wide range in slopes. 

The mean concentration C for the entire flow and Cy, the mean concen- 
tration in the region below dy, have been plotted in Figs. 6(c) and 6(b) in 
terms of the discharge and for each slope. It is interesting to note that the 
mean concentration for the entire flow is nearly constant for a quite wide 
range of discharges but increases very considerably with the slope. The 
mean concentration Cy for the region below dy has similar characteristics. 
In Fig. 6(a) the concentration C-y at the transition depth is plotted. This 


concentration varies but little with the depth but increases with the channel 
slope. 


Analysis of the Data 


Distzibution of Air in Aerated Flow 


A qualitative concept of air-entrainment phenomena has been established 
with the aid of the air-concentration determinations made transverse to the 
direction of self-aerated flows, high-speed photographs and motion pictures 
of the flow, and other similar information. Since it has been shown that aer- 
ation of flows on a steep slope begins in a region where the boundary layer 
generated by the channel bed has reached the water surface, it is to be pre- 
sumed that aeration depends primarily upon the turbulence intensity of the 
flow. Aeration appears to occur when the transverse velocities of turbulence 
are sufficiently strong near the air-water interface to cause clumps of water 
to break through the surface into the air and then fall back by gravity into the 
flowing stream. This breaking away of clumps of water from the main stream 
and the falling back into the stream with attendant splashing and breaking into 
a heterogeneous spray of globules and droplets is associated with the insuf- 
flation of air into the stream. It occurs particularly with water flow at high 
velocities on steep slopes. The air carried back is then distributed through- 
out the flow by turbulent transfer. Observations by means of instruments indi- 
cate that two regions of self-aerated flow develop: (1) an upper region of 
heterogeneous clumps, globules, and droplets of water ejected from the flow- 
ing liquid stream into the atmosphere at more or less arbitrary velocities, 
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and (2) a lower region consisting of air bubbles distributed through the flow 

by turbulent transport fluctuations, which can be described by some boundary- 
layer equation. Between the two regions is a transition zone defined by a 
transition depth which is a fluctuating surface necessarily at a statistical 
mean elevation above the channel bottom. 


Distribution of Air in Upper Region 

If it is assumed that the components normal to the bed of the turbulent 
velocity fluctuations are randomly distributed such that the mean value is 
zero, those in the outward direction will constitute one-half of a Gaussian y y 
distribution. Those in or near the surface that penetrate the surface will 3 | 
carry large clumps and smaller globules of water into the atmosphere against 4 
the force of gravity a distance proportional to the square of the individual 
transverse velocities. The maximum height of the trajectories that they 
follow will also be randomly distributed about zero. (The path will be a long, 
flat trajectory because of the large, longitudinal component of the local ve- 
losity.) If the distances from the surface from which the particles are pro- 
jected are represented by one half of a Gaussian distribution, the frequency 


f(y') of particles projected a distance y' above the transition depth may be ex- 
pressed as 


2 h 

(4) 

h./r 
where h is a measure of the mean distance the particles are projected above 


dy. The factor 2 is used to indicate that only particles being projected in the 
outward direction are being considered; that is, half of the complete Gaussian 
co 


is being taken so that J f(y') dy' = 1. Furthermore, since it is presumed 


that the flow is uniform on the average, particles are equally likely to be pro- a 


f(y') = 


jected at all points of the transitional surface. Then the proportion of all 4 
particles leaving a unit area of the transitional surface that reach or pass : 
through a corresponding area at elevation y' will be is: 


2 
h 


and the number of particles reaching or passing through this area in both di- 
rections during the sampling period will be 


Nv = N d 


(6) 


where Ny is the total number of particles that leave or return to the transition 
surface per unit area during the sampling period. The average number of 2 
particles of water that reach or pass through any horizontal area per unit ie 
time may be assumed proportional to the average concentration of water parti- 
cles. Therefore, using Eqs. (5) and (6), the air concentration at any distance 
y' above the transition level in terms of the concentration at the transition 
level Cy and the representative height h is 


= 
— 
| 
a 
i 
on 
Py ad y (5) 
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(7) 
dy 


I- Cy 


Eq. (7) is, of course, the cumulative Gaussian probability curve and is appli- 
cable to the aeration phenomenon only for y' 70. The gradient of the concen- 
tration is from Eq. (7) 


- 


dy 


de. 
and when y' = 0, —— is a maximum whose value is 


dy 
dy /mox (9) 


and occurs where the concentration is Cp and the depth is dr. 

The constants C-y and h are properties of the flow related to the intensity 
of the turbulence in the neighborhood of dy which, in turn, is probably de- 
pendent upon the velocity and roughness of the channel. 


Distribution of Air in Lower Region 

The air entrained at the transition surface by the turbulent fluctuations is 
then distributed throughout the flow in the region below dy by the turbulence 
in the stream resulting in a statistical equilibrium between the bouyancy of 


the air and the concentration gradient. Such equilibrium is described by the 
well-known equation 


CV + (10) 
where C is again the air concentration at a normal distance y from the bed; 
Vp, the rising velocity of the air bubbles, is taken as being negative; and €) 
is a mixing parameter for the air-bubble transfer. In aerated flows the value 
of €) as a function of y is unknown, and a number of assumptions as to the 
form of the function can be made. However, it may be presumed that the air 
bubbles are transported in the same manner as momentum; and it appears 
reasonable as an approximation to assume €p to be proportional to €,,, the 
momentum mixing parameter, as is done in the case of transport of other 
substances such as sediment. In open channels in the nonaerated condition 
where the shear is linearly distributed from zero at the surface to a maxi- 
mum at the bed and the velocity profile is considered to be logarithmic 
throughout the depth of flow, it has been shown(9) that € m has a parabolic 


form being zero at the surface and at the bed with a maximum at the center- 
line and can be expressed in the form 


moe 
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where VT,/p = V,, the shear velocity, and 7, is the boundary shear at the 
bed, p is the density of the fluid, Bisa proportionality factor in €, = Bép’ 

and k is the Karman universal constant. In the case of transport in a closed 
channel, however, the shear at the upper boundary is not zero and the maxi- 
mum velocity occurs at a point between the upper and lower boundary, at 
which point the shear is also zero. This requires that €,, also be zero at this 
point and that the distribution of €,, consists of two parabolas having zero at 
the bed, the point where the maximum velocity occurs, and at the upper 
boundary and maximum values at intermediate points above and below the 
point of maximum velocity. It has been shown(10) however, that suspended 
matter is actually transported transversely across the level of zero mixing 
coefficient, and it has been suggested that the mixing coefficient be considered 
constant in this region. Since it is likely that the vertical transport of sus- 
pended matter depends only upon the transverse velocity fluctuations in the 
vertical direction, such transport may well occur even if the correlation be- 
tween the vertical and longitudinal fluctuations approaches zero. 

Measurements of velocity distribution in aerated flows show(4,5) that the 
velocity tends to decrease as the upper limit of the flow is approached and 
that the maximum occurs well down in the flow proper. This implies that a 
shear is developed in the transition zone. A plausible explanation for such 
a shear is the change in momentum engendered by the return of water clumps 
and droplets to the flow at the end of their trajectory through the atmosphere 
above dy. In the course of that trajectory, they suffer a loss of velocity due 
to atmospheric drag on the mass. The shear developed at dy is analogous to 
an upper boundary and a distribution of the type found for closed channels 
might be expected. 

The complex distribution found for closed channels may be approximated 
by a parabolic distribution over the entire region below dy, with the maximum 
value incorporated in the factor 8. This approximation is equivalent to as- 
suming a distribution for €p as given by Eq. (11). If this value for <p, given 
in Eq. (11), is substituted, Eq. (10) becomes 


which upon integration, assuming that Vp) is independent of y, is 


= y 
( (14) 


in which 


Vb 


Bk Vy (15) 


d 
and C, is a constant whose value is the concentration at y = — 
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Application to Experimental Data 


The use of Eqs. (7) and (14) as a description of the air-concentration 
distribution characteristics depends upon the degree to which they fit the ex- 
perimental data. There are two factors which have a bearing on this: the 
form of the equation and the magnitude of the constants. Whether or not the 
form of the equations is reasonably correct depends upon the extent that the 
assumptions made in their development approximate the actual mechanism. 
The constants involved depend upon the intensity of the turbulence and at the 
present cannot be determined analytically but must be evaluated empirically 
from systematic experimentation. 

Eq. (9) permits the determination of the transition depth dy and the concen- 
tration Cy at that depth. From a plot of the concentrations in terms of the 
normal distance from the bed, the point where the gradient = is a maximum 
was graphically located, as shown in Fig. 7(a). The corresponding values of 


dp and Cy are noted. Using this value of Cy, the ratio rte was plotted 


as a function of y on so-called probability paper, upon which a cumulative 
Gaussian distribution plots as a straight line. The value of h, or rather the 
value of o, the standard deviation of the distribution, wherein h = oy?, is ob- 
tained from the slope of this straight line. Alternately, h can be determined 
directly from Eq. (9) by evaluating the gradient at the point where - isa 
maximum and inserting the value of Cy obtained for that point. Fig. 7(b) 
shows a plot of the concentrations in the upper region of the flow plotted in 
this manner with the standard deviation indicated. It is apparent from this 
plot that the concentrations possess a cumulative Gaussian distribution. 

In the region of flow below dy, Eq. (14) indicated that the concentration C 


is a function of ec which, when plotted logarithmically, results in a 


straight line whose slope is equal to the exponent z and C, is the value of C 

at Pa = 1. The concentration data for the experiment have been plotted in 
this manner in Fig. 7 to show the degree of agreement with the form of Eq. 
(14) and to indicate the evaluation of C, and z. Inasmuch as C approaches 
infinity as y approaches dy, the experimental data must depart from the curve 


for the larger values of ee but in this case the data agree with the curve 
up to = 0.9. 
dy 


The straight lines drawn through the data in Figs. 7(b) and 7(c) have been 
transposed to the original plot in Fig. 7(a) to show the characteristics of the 
curves on the arithmetic scales. Eqs. (7) and (14) have also been plotted 
along with the data for a number of typical experiments using various dis- 
charges and slopes shown in Figs. 8 and 9 for purposes of comparison. Be- 
cause o*.the character of Eq. (14) in the neighborhood of y = dy, the two 
curves are connected by a transition to pass from one curve to the other. 


Relation of Parameters to Flow Conditions 


Based upon the above analysis of the aeration phenomenon, a number of 
parameters were developed as being descriptive of both the magnitude and the 
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distribution of the air. These parameters varied with the flow conditions 

such as the channel slope and discharge, which in turn, along with the rough- 
ness characteristics of the channel, govern intensity and scale of the generat- 
ed turbulence. The flow turbulence is created initially at the bed by the wakes 
and eddies formed by the flow over the roughness elements; the turbulent 
eddies are then diffused upward into the flow stream. The intensity of turbu- 
lence at the transitional surface which appears to be paramount in the aeration 
process should then depend upon both the initial generation and the depth of 
the mixture to the transitional surface. The intensity of the turbulence is re- 
lated to the boundary shear T,, on the bed, which can be expressed by a so- 
called shear velocity Vep = V7T>/p. Because of the shear that is presumed to 
exist at the transitional surface, it is not possible in the present experiments 
to differentiate the bed shear from the total shear. As a first approximation, 
the total shear or a shear velocity based upon the total boundary shear may 
be used as a measure of turbulence intensity. To test this hypothesis, the 
mean concentration was plotted in terms of V,/ay2/3 where V, is defined as 


= (16) 


and dy is the transitional depth previously defined. The expression V,/az?/3 
is empirical and was chosen as that form which best correlated the data. 
This plot is shown in Fig. 10. The measured values of the mean concentration 
fall reasonably well on a single curve with a tendency for the points corre- 
sponding to the smaller discharges to depart from the mean curve. When the 
mean concentration in the area below dy, C7, is considered as being the air 
which is insufflated and distributed through the flow, the correlation is even 
better, as shown in Fig. 11. For the higher discharges for each slope the 
points all fall on a smooth curve, while there is a clear departure from the 
curve for the low discharges on the corresponding slopes. It appears that as 
the discharge decreases below a certain value for a given slope, an additional 
factor becomes important in reducing the concentration. The nature of this 
factor is not readily apparent, but it should be noted that a decrease in dis- 
charge corresponds to a decrease in dy, and it is quite possible that an insta- 
bility in flow pattern arises when the depth becomes small. Considering the 
data as a whole, it is rather significant that when plotted in this manner, mean 
air concentrations ranging from about 0.05 to 0.75 corresponding to discharges 
ranging from about 4 to 15 cfs on slopes with ranges from 7.5 to 75 degrees 
are correlated along a single curve. 

The parameter V,/dy2 3 includes the depth dy to the transitional surface. 
This depth is a function of the air-concentration distribution and is difficult 
to determine explicitly. It is, however, related to the discharge and slope and 
in general to the channel roughness. In these experiments the channel rough- 
ness is a constant so that for these data, dy depends only on the slope and 
discharge. In terms of these measured characteristics, v,/az2/3 is related 
to S/q1/5 where S = sina and a is the angle of inclination of the bed and q is 
the unit discharge of water. Consequently, the air concentration that may be 
expected in a high-velocity open channel can be related to the flow charac- 
teristics. This relationship is shown in Fig. 12, where the mean concentration 
as measured is plotted as a function of S/q1/5 for the particular channel in 
which the experiments were performed. For channels of different roughness 
characteristics, the curve will probably be similar but shifted an amount to 
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correspond to the different roughness. A similar plot is given in Fig. 13 
showing the mean concentration between the bed and the transitional surface 
where the depth is dy. Here again, as in Fig. 11, correlation between Cy 
and S/q1/5 appears to be better than that between C and S/q1/5. 

The results shown in Fig. 13 provide a means of estimating the mean air 
concentration that may be expected in an equally rough, wide open channel of 
known unit discharge and slope. The depth of flow or the magnitude of the 
various previously described depth parameters, which are important in de- 
sign, is a function of the air concentration. 


Effect of Entrained Air on the Flow 


It was suggested above that the air insufflated into the flow is distributed 
throughout the flow by turbulent transport. The entrained air has the effect 
of changing the flow from one of water alone to the flow of a mixture of air 
and water. Several characteristic depths were defined, such as d,, to repre- 
sent the uppermost extent of the water particles projected from the flow, d 

to represent the depth of flow of the water alone, and dy to represent a new 
equilibrium depth brought about by the presence of air in the flow. For de- 
sign purposes, it is of considerable importance to know how these parameters 
and the actual flow velocity are related to the bulk flow properties or what 
differences might be expected in these parameters in relation to a corre- 
sponding nonaerated flow for which the depth may be computed by means of 
well-established flow formulas such as the Manning or Chezy formula. 

In order to investigate such a relationship, experiments were performed 
in the experimental channel to establish the resistance of the channel to non- 
aerated flow, that is, flow of water alone. The slope was decreased so that 
for comparable discharges the flow would be nonaerated. These slopes 
ranged from approximately 1 degree to a maximum of about 6.5 degrees, 
above which aeration was observed to occur. The discharges were in the 
same range as those used in the aeration experiments. Normal flow was es- 
tablished for each combination of slope and discharge, and the mean depth 
was determined by averaging the depths measured at various points along the 
centerline. 

The relationship between the depth of flow for the nonaerated condition and 
the discharge and slope is shown in Fig. 14, where the measured depth d,, is 
plotted as a function of q/s1/2 on logarithmic scales. The equation of the 
line drawn through the points is 


90 (sina)! (17) 


For the range of the variables encountered in the nonaerated flow, it appears 
that the Chezy formula best describes the data and the coefficient 90.5 charac- 
terizes this particular roughness. The mean velocity is then 


Vm= 


If it were assumed that the same type of flow equation applied to aerated 
flow, the mean velocity of the air-water mixture flowing at a depth dy is 
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V = “(sina)”? (9) 


where V is the mean velocity of the air-water mixture and C is the Chezy 
coefficient. Assuming that the velocity of the air and water components are 
equal, then Eq. (19) can be written in terms of the water discharge, as 


= Cdd, (sina) (20) 


where d is the mean depth computed by Eq. (1) and represents the cross- a 
sectional area of the water. Eq. (20) was tested by plotting the experimental 
data for aerated flow in the form ddy 1/2 as a function of q/(sina)!/2 in Fig. 
15. Superimposed upon the plot is the line for nonaerated flow in the experi- 
metal channel taken from Fig. 14. The data cluster about the straight line 

for nonaerated flow, indicating that an equation similar to Eq. (20) will corre- 
late the data, and, in addition, when the parameters d and dy are used, the 
Chezy resistance coefficient is essentially the same as that for nonaerated ie: 
flow. Fig. 15 shows that at least for the conditions of these experiments, air- ie 
entrained flow, when considered to be a flow of a mixture, follows the same 
laws as for nonaerated flow when depth parameters related to the properties dl Mi 
of the mixture are used to define the linear flow dimensions. The agreement \ 
of the data for aerated and nonaerated flow shown in Fig. 15 substantiate 
somewhat the deductions made in regard to the air-entrainment process, 

An evaluation of d, dy, and dy for specific flow conditions is not possible 
from Fig. 15. A clearer idea of the changes in these parameters will be ob- 
tained if they are compared with the depth d,,, computed in the usual manner, 
using Eq. (17) which includes for this channel the Chezy coefficient of 90.5. 

It is the depth at which the water would flow in this channel at the given dis- 
charge and slope assuming no air entrainment. The ratios dy/dm, dp/dm, 
and d/dy as functions of mean air concentration are plotted in Fig. 16(a), 
16(b), and 16(c). In all cases, of course, the ratios approach unity as the 
mean _concentration approaches zero. That is, for water flow only, 

d,, = d = dy =dm. As the mean air concentration increases, the ratios depart 
from unity in such a way that as C approaches, 1.0, d becomes very small, 
while dy and dy necessarily increase. 

The upper limit of the spray d, increases very considerably with increased 
air concentration which in turn means an increase in the turbulence intensity 
to project water particles farther into the atmosphere. For the highest 
concentrations measured, the value of d, is nearly four times as large as the 
nonaerated depth for a corresponding discharge and slope. The transition 
depth, as might be expected, also increases relative to the computed depth 
since increasing concentration bulks the flow and makes the effective depth Bik 
greater. On the other hand, the mean depth decreases with increasing concen- 
tration. It appears from Fig. 16 that concentration has relatively little effect 
on dy and d until the air concentrations are greater than about 50 per cent. 
For the slopes where the air concentrations were the greatest (being about 85 
per cent), d was reduced to about one half of d,, while dp was nearly twice 
as large as d,, 

The decrease in d with increase in mean air concentration implies a corre- “s 
sponding increase in the mean water velocity above that computed on the basis 4 
of dy, since d is the cross-sectional area of the water flow. The trend of the q 
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mean water velocity in terms of the nonaerated velocity in this channel is 

shown in Fig. 17. Here again, the velocity of aerated flow is about the same 
as that of a corresponding nonaerated flow until the mean concentration is in 
the neighborhood of 50 per cent. It then increases rather rapidly. 


CONCLUSIONS 


1. Self-aerated, high-velocity, open-channel flow appears to consist of 
two distinct phases: an upper region consisting primarily of independent 
droplets and larger agglomerations of water that move independently of the 
stream proper, and a lower region in which discrete air bubbles are suspend- 
ed in a turbulent stream and are distributed by the mechanism of turbulence. 

2. In the upper region, the distribution of water-air agglomerate and drop- 
lets can be adequately described by the Gaussian cumulative probability 
equation, the constants of which are properties of the flow conditions. In the 
lower region, the distribution of the air agrees closely with an equation for 
turbulent mixing based upon an approximation for the distribution of the mix- 
ing parameter. 

3. The magnitude of the entrained air concentrations apparently depends 
upon the intensity of turbulent fluctuations generated at the rigid bed and the 
depth of flow. The mean concentration could be correlated with the parameter 
V,/ay2/3 or in terms of the flow characteristics; the mean concentration C 
is a function of S/qi/5, where S is the sine function of the slope angle, and q 
is the unit discharge. 

4. The maximum height of water (or spray) above the rigid bed relative to 
a corresponding nonaerated flow increases rapidly with mean air concen- 
tration in the flow. The transition depth is greater than the depth of a corre- 
sponding nonaerated flow because of the bulking effect of the entrained air. 

5. The mean velocity of an air-entrained flow is greater than that of a 
corresponding nonaerated flow by an amount that increases with the air 
concentration and corresponds to a decrease in the mean depth. 

6. When the effective depth of aerated flow is taken to be the transition 
depth and the mean depth is the cross-sectional area per unit width of channel, 
the flow is described by a formula of the Chezy type involving approximately 


the same roughness coefficient as applies for nonaerated flow in the same 
channel. 
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LIST OF SYMBOLS 


- Concentration, ratio of volume of air to volume of air plus water. 


ara 


- Mean concentration in vertical section. 
C, ~- Concentration at point y = d7/2. 
Cr - Mean air concentration in region between lower boundary and dy. 


Cy - Concentration at transition depth y = dy. 


d - A mean depth of flow defined by Eq. (1). 

di. - Upper limit of flow or value of y where C = 0.99. 

d; - Transition depth between upper and lower regions of flow. 
din - Mean depth in nonaerated flow. 

g - Acceleration due to gravity. 
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h - A mean height to which water particles are projected above dp. 
k - Von Karman universal constant for velocity distribution. 
Ny: - Number of particles reaching or passing through a horizontal area 


a distance y' above transition depth dy. 


Ny - Number of particles that leave or return to an area at the transition 
depth. 
Py: - Proportion of all particles leaving area at transition depth that 


reach or pass through horizontal area at y' above transition depth. 


- Total-water discharge. 


- Water discharge per unit width of channel. 


- Mean velocity in vertical section = q/d(d). 


Shear velocity = g(d7) sina. 


Vp ~~ Rising velocity of air bubbles. 


Q 
q 
V - Local velocity. 
Vv 
Vv 


- Normal distance from channel bottom. 


- Outward normal distance above transition depth. 


Zz - Exponent in air-concentration equation. 

a - Angle of inclination of channel. 

B Constant. 
bb - Mixing coefficient for air bubbles in turbulent flow. ‘ 
San - Mixing coefficient for momentum transfer. 


- Kinematic viscosity of water. 


- Density of water. 


Vv 
p 
o - Standard deviation of air distribution above dy. 
T = Local shear force. 


- Boundary shear force. 
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Fig. | - Variable Slope Channel for 
Air-Entrainment Experiments 


Fig. 3 - Probes for Air-Concentration 
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Fig. 2 - Artificial Roughness Elements Installed on 
Channel Bed 
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Fig. 4 - Comparison of Air-Concentration Profiles 
at Stations 35 ft and 45 ft From Inlet 
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a 
d 
‘ 
ASCE 
Mee 
9 
am 
a 
ay | | 
8 40 
= 
| 
~ 
q 
~o. 
> 
— 
a 
a ~ 
a 7 
: 


1890-24 


Concentration , C 


Mean 


HY 7 December, 1958 


Symbol Slope Degree 
75 
15 
225 


S/q* 
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0.565 3uel 0.113 040512 14502 | 
0.550 350k 1.852 0.113 0.068 1.427 | 
0.555 37.7 1.977 0.087 0.078 1.423 | 
0.555 40.0 2.16 2.092 0.079 1.389 | 
0.525 06 2.2 0.113 0.062 14372 | fe 
0505 0.490 2034 0.122 0.080 1.289 
0.453 400 2643 04090 1.218 | 
| Whe 3.398 1.8 2043 0.151 0.088 1.208 
e352 0.326 42.0 2.450 0.234 04095 16157 | 
78 0.613 37.0 1.539 0.107 0.053 1.737 Be 
640 £0 0.630 1.673 0.102 0.061 1565 | 
M654 0680 0.091 0.061 1.636 
0.82 0.645 45.7 2430 0.098 0.073 1.606 | 
0.81 0.6657 48.2 0.0806 0.090 14570 | 
0.60 0.640 2458 0.086 9.095 
0.78 0.610 L7.0 2.69 0.087 0.0968 
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